ASYMPTOTIC DISTRIBUTION OF THE ROOTS OF THE EHRHART 
POLYNOMIAL OF THE CROSS-POLYTOPE 



MIGUEL RODRIGUEZ 

Abstract. We use the method of steepest descents to study the root distribution of the Ehrhart 
polynomial of the d-dimensional cross-polytope, namely Cd, as d — >■ oo. We prove that the 
distribution function of the roots, approximately, as d grows, by variation of argument of the 
generating function X]m>o C.d{m)t"^'^^~^ = {1 + 1)'^ {1 — t)~'^~^t^~^ , as t varies appropriately on 
the segment of the imaginary line contained inside the unit disk. 



1. Introduction 

There has been of recent some interest in the root distribution of Ehrhart polynomials. In [5], 
Beck et al. come up with bounds for the coefficients and roots of Ehrhart polynomials. They 
found them to be contained inside a disk of radius d\ centered at —1/2, where d is the degree 
of the Ehrhart polynomial in question. This result was subsequently improved upon by Braun 
in [4j and [5J, and Bey et al. in [3], where the radius of the above mentioned disk is reduced 
to the order of d'^. All of these last results depend Stanley's non- negativity theorem, which 
states that the coefficients of the Ehrhart polynomial are, when the polynomial is expressed in 
a binomial basis, positive. In [7J, Vaaler et al. prove that the roots of the Ehrhart polynomial of 
the cross-polytope all have real part —1/2. In [16j, Rodriguez- Villegas generalizes the last result 
by proving that if the generating function Y2rri>o'^ i^)^^ ^ where V is some polynomial, has all 
of its roots on the unit circle, then, V has all of its roots on a vertical line in the complex plane. 
These last results depend, not on the positivity of the coefficients, but on functional equations 
that are satisfied by the above generating function. Here, we look at the zero distribution of 
the Erhart polynomial of the d-th dimensional cross-polytope, namely C^, as d grows. The 
form of the generating function allows us the use of the method of steepest descents to study 
the asymptotics of Cd on the line {z £ C : Re{z} = —1/2}. Since our aim is to obtain an 
asymptotic formula which approximates Cd uniformly on ever large segments of a vertical line, 
we will have to deal with eventualities that arise in the use of the method of steepest descent for 
uniform approximations; we will have to deal with coalescing saddle points and with a saddle 
point coalescing with a singularity of the integrand. The terminology will all be explained below. 

The d- dimensional cross-polytpe is the convex hull in M'^ of {ibej|j = 1, . . . ,d} where ej = 

{Sji, . . . , 6jd) where 6ji is the kronecker delta. Let Cd be the function Cd{m) = n mVd) 

for m £ 7j-^. Cd is actually a polynomial, called the Ehrhart polynomial of the d- dimensional 

cross-polytope. For introductory material on the cross-polytope and Ehrhart polynomials see 
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|. Cd can be expressed as a sum of binomials: 

d 

(1.1) Cd{x) = ^ 



3=0 



d 



j + x 
d 



The sequence {Cci{m)}rn>o has generating function: 
(1.2) 8d{t) 



m>0 



(1 - t)rf+l 



We can also write: 
(1.3) 



1 

2TTi 



£d{t)t 



-x-1 



dt, 



where F C C is some simple closed curve, oriented counterclockwise with respect to the bounded 
connected component of its complement. F contains 1 in the bounded connected component of 
its complement, while is in the unbounded connected component. We can derive the above 
formula as follows. 

{i + tY 



£d{t)f 



-x-1 



dt 



(1 - tY+^ 



r'^-^dt 



d 

E 

j=0 
d 

E 

j=0 



d 
j 



-x-1 



-dt 



27riRes 



j=0 



-1) 



d 



d+l 



X 
d 



j + x 
d 



= 2mCd{x) 

The above integral representation allows us the use of the method of steepest descents to produce 
an asymptotic formula for Cd- An explanation of the method of steepest descents and our use of it 
will be given in the sections below. In these notes we use Landau's asymptotic notation. If / and 
g are complex valued functions with domain C C, we write f{z) = 0{g[z)), or f{z) <^ giz), 
as z — 7- a, if there is some constant C > such that |/(z)| < C|5((2;)| for all z sufficiently close 
to a in D (when a is oo, by all z sufficiently close to a we mean for all sufficiently large |2:|). 
The previous situation can also be described by the symbols g{z) = Q{f{z)) or g{z) » f{z) 
as z — )• a. Sometimes we will omit "z — )• o" if we believe that its presence is obvious from the 
context. By f{z) = o{g{z)) as z — )• a we mean — )• as z — )• a. By f{z) = uj{g{z)) as z — )• a 

we mean — >■ +oo as z — t- a. By f{z) ~ g{z) as z — t- oo we understand ^||y — t- 1 as z — t- a. 
We will use f{z) = Q[g{z)] to denote that g{z) <^ f{z) <^ g{z) as z — a. We will use D to 
denote the closed unit disk. Our asymptotic formula for Cd is given in the following theorem. 
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Theorem 1.1. Let e £ (0,1/2). Let Sd be the sequence of sets Sd = {z £ C : Re{z} E 
[e, 1 — e], |2;| < d — Vd}- Given x G Sd, define a = a{x) = —i (^ i^^y — i^^y^ ~ • Cm/fn 
be the m-th degree coefficient of the Taylor series o/ log{iSd(t)(l — t)f^^^^} around a. Then, 
there is a complex valued function T{d, x) defined for x £ Sd, such that, for large enough d: 



(1.4) Cd{-x)=T{d,x) + {-l)'^+^T{d,l 

and, 



^[(2(i + l)-^r(x) + 0[(2d + l)-^-i] ; if x = 0{.J\ogd) 



(1.5) F{d,x) , 

g£:rf(a)a^Xd(x){l + 0{lm{xY^'^^)} ; if ^\x\<.d-^ 
Where D is some positive constant, and. 



Ux) = |^^^exp|i^^F(T)|c/'(T)dT 



1-1/3 

,u,here FiT) = 2T-^^^^^, UiT) = T(^,+^). 

Proof. The roles that F{T), U{T) play in our use of steepest descent, and, in our study of the 
roots of Cd will become clear in the proofs, which will all be given below. □ 

The above result might be stated briefly as: "The function Td is approximately given by the 
generating series £d{t)t^ evaluated at the saddle point a;(x)". We will be able to determine the 
root distribution of Cd on the line {z £ C : Re{z} = —1/2} just by measuring the variation of 
the argument £dit)t^ , for t varying in the intersection of the imaginary axis and D. Theorem 
|1.1| only refers to the asymptotic nature of Cd in a restricted range which, as we will see below, 
contains all the roots. 



2. Asymptotic distribution of roots of Cd 



2.1. Roots with imaginary part in [—^/logd, \/log d] . We first study the roots of Cd , which 
have imaginary part in [— -v/log d, \/log d]. For large enough d, and for x = 1/2 + ir with 



r £ [0, \/log d] , we can write, using (4.11): 



vr f 2Re{(2d+ l)-^r(x)5d(a;)} , if d is even 
Cdi-x) = < 

smvrx \2ilm{{2d+l)"-'T{x)gd{x)} , if d is odd 

where gd{x) = 1 + (k^) • Notice that gd is analytic in Sd = {x : Re{x} £ [e,l — e],Im{x} £ 
[— Vlog (i, 0]}. We already know the roots of Cd have real part —1/2. Hence, to find roots of 
Cdi—x), we solve: 



(2.1) 



arg{(2d + l)-^T(x)gd(x)} 



^ mod vr 
mod IT 



if d is even 
if d is odd 
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Then, using Stirling's formula, as c? — )• oo, 
arg{(2d+l)-"r(x)5d(x)} 



arg{(2d+l)-^r(x)} +arg[5rf(x)] 
^logTTTT^^^ +0(1) 



e{2d+l) 

Here we take the principal branch of the argument function; that for which arg(l) = 0. 

Now, if / = n + it; is some analytic function with real and imaginary parts n, v, we can write 
the partial derivatives of the component functions u and v in terms of the derivative of /. In 
fact, if / is a function of the complex variable x = cr + ir, with a and r in M, then, using the 
Cauchy-Riemann equations for example, 

d 



dr 



v{x) = Re{nx)} 



This means that we can use Cauchy's theorem to bound the partial derivatives of u and v in 
terms of /. Then, since gd is analytic and different from inside 5^, we can use Cauchy's 
estimate (see [I]) to obtain a uniform bound for the derivate of the for all d. So, for large 
enough d, 

d 

dr ' Tr^'^^^ 



^arg {(2d + l)~-T{x)ga{x)] = ^{2d + l)--T(x) 



dr 



rlog- 



+ 0(1) 



e{2d+l) 
= log(2d+l) + 0(l) <0, 

In the estimates above we used that < logd. We also took the argument to be the 
principal branch with arg = on the positive real line. It follows that for large enough d, 
arg{(2(i + l)~^r(a;)(7rf(x)} is a monotone function of r for x G 5 (this method is used by 
Lagarias to count the zeros of differenced L functions in [E]). Let a and b be such that 
— log d < a < b < log d. Then, by (2.1 ) and the work done above, as d — )• oo. 



(2.2) #{x £ Sd ■ X is root of £^ with Im{x} € [a, b]} 
1 



TT 



arg <^ {2d + 1) 



1 



+ ib] }+Os{l) 



Where the error above is actually less than l.Then, as d - 
#{x £ Sd ■ X is root of Cd with Im{x} € [a, b]} 



+oo: 
logd 



(6 



, as d 



oo 



2.2. Roots with absolute value in [^/logd,d — \^]. We now study the roots of Cd{—x) 



which have imaginary part in [-v/logd, d — ^fd\ . For large enough d, and for x 
T G [Vlog d, d - 



I + ir with 



we can write, using (|l.4|): 



2ilm| |^f(i(a)a^X^(x)5(i(x)| , if d is even 
2Re{gfrf(Q)a^Xd(x)5rf(x)} , if d is odd 
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where Qdix) = 1 + 0(r^^/^^). As in the previous section is analytic in Sd = {x : Ke{x} E 
[e, 1 — e],Im{3;} G [—d+ v^, 0]}. Also, as in the previous section, to find roots of x), we 
solve: 

f 3C9 1 TT mod vr , if (i is even 

(2.3) arg -^£,(a)a^I,{x)gd{x) \ = {^ 

I ^'^3 J mod vr , if d is odd 



Now, U'{T) is in the first quadrant. In fact, looking at (4.25), it is not hard to see that 



Im{C/(T)} = i, while Re{U{T)} > Then (see (4.15) for the definition of Km), 



n/3exp{-|r§F(r)}dr ^ 

(2.4) arg[X,(x)] = arctan — ' ^ ^ G 0, - 

n/3exp{-|r^F(T)}Re{C/'(T)}dr ^ 3 

Then, as d — )• 00, 

f 3C 

arg I ^^£d{a)a''Id{x)gd{x) 
3C 

= argTTTT + + 1) arg(l + q) + arg(a^) + arg[Id{x)] + aTg[gd{x)] 
= -{2d + l)arctan(|a|) + rlog|a| +0(1) 
Use Cauchy's estimate to obtain, for large enough d, 

|:«g{||f.(a)»^X.(.to(.)} = - ^} +0(1) 

It follows that for large enough d, avg ^'^i^Sd{a)a^Id{x)gd{x)^ is a monotone function of r. 
Now, notice that when \x\ = 0(-y/log d), then, arg | |^<5d(a)a^Xd(x)| — arg{(2(i + l)~^r(x)} = 
0(1). Then, if-d+ ^ < a < b < d — \fd^ as d — >• +00: 

1 

#{x G S'rf : X is root of £d with Im{x} G [a, 6]} ~ — {arg[J^(d, 1/2 — i6)] — arg[J^((i, 1/2 — ia)]}, 

vr 

as d — )• 00. 

2.3. The root of largest absolute value. Again, 
f 3C 

arg I -^£d{oL)a''Xd{x)gd{x) 
3C 

= argTTTT + (2d + l)arg(l + a) + arg(a''') + arg[Xrf(x)] + aTg\gd{x)\ 
= -| - (2d+ l)arctan(|a|) +rlog|a| - ^ + arg[Xd(x)] + 0(r"^/2^) 



Assume d is odd. The d even case is similar. Equation (2.3) allows us to estimate all the roots of 



Ld that are of absolute value 0(d— \^), the range of validity of our approximation. Assuming 




Figure 2. Error from the counting function of the roots of £500 on [0, 500]; i.e. 
#{x G 5500 : X is root of Cr,oo with Im{a;} G [0,r]} - ^{arg[J"(500, 1/2 - ir)] - 
arg[J"(500, 1/2 - iO)]} for r G [0, 500]. 



that the largest root of Cd[—x) is x 
solving arg||g^£:rf(a)Q!^Trf(x)5rd(a;)| 



1/2 + ir = 0{d — v^), we could obtain an estimate by 
-[(d-l)/2]7r, or, 



(2.5) 



{2d + 1) arctan(|a|) + r log \a\ 



2d -5 



7r-arg[Xrf(x)] + 0(T-V28) 



In fact, Cd is a polynomial with real coefficients and d distinct roots. This last bei ng a conse- 
quence of the monotony of arg ||^£^(;(a)a^'X(;;(x)(7d(x)|. So if we manage to solve (2.5) in the 
desired range, we will have found the largest root. By ( 4.13| ), if r = o{d) as d ^ 00, then, the 
RHS of the above equation is o{d). So it must be that r = Q{d). Assume, for the moment, that 
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T = d — f{d), where f{d) = and positive. Let y = f{d)/d. Then, using (4.13), we may 



write (2.5) in terms of y: 

1 16d+l 



^y + ^^^y' + 0{dy') = -vr - arg[X,(x)] + 0{d-'/'') 



By (2.4), ^TT — arg[Xrf(x)] is positive and bounded away from zero. It then fohows that: 



,or, 

/n\ 1/3 f o ^ 2/3 

fid) = ^ (gj I 2^ - arg[X,(x)] + 0(d-i/2«)| 
In any case we have shown that the largest root has imaginary part r, such that d — T = Q{\/d). 

Table 1 . r is the imaginary part of the largest root of Cd 



d 


r 




d-T 


100 


91.9987057014 


1 


7238266002 


200 


189.7372321215 


1 


7549086218 


300 


288.1562327578 


1 


7692238245 


400 


386.8992027271 


1 


7780517454 


500 


485.8385218444 


1 


7842344425 


600 


584.9118679958 


1 


7888958567 


700 


684.0835726177 


1 


7925842603 


800 


783.3310874715 


1 


7956041698 


900 


882.6391445854 


1 


7981404758 


1000 


981.9968699646 


1 


8003130035 



3. Proof of the asymptotic formula. 

Let {xd}d>o be a sequence of complex number such that Xd = (7d + ird, where ad, Td G M. Let 
e G (0, 1). Assume ad G [e, 1 — e] and \Td\ G [0,d — \^]. In this section we use the method of 
steepest descents to study the asymptotic value of: 



(3.1) J^£d{t)t-'=''-^dt 



as d — 7- +00, where F is as in (1.3). The basic result needed for the implementation of the 
method is Watson's lemma, which we subsequently introduce. Most of the method of steepest 
descent consists in modifying integration contours in a way that allows for the use of the lemma. 
From now on we drop the subindex in most of the quantities that vary with d. We hope this 
last will make our notation clearer. 
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4. Watson's lemma 

A basic result needed for the application of the method of steepest descents is Watson's 
Lemma (see [15]), which deals with integrals of the form: 

(4.1) re-^'^Wg(t)dt 



as X —7- +c«, where (j){t) is some complex valued function of the form t + O(t^) as t — t- 0"*" that 
has an increasing real part on [0,a], and, q{t) ~ as t — )• O"*". A typical version is: 

Theorem 4.1. Let q{t) be a function of a positive real variable t, such that 

g(t) ~ ^ a,t("+^-'^)/^ {t^O) 

s>0 

where A and fi are positive constants. 

(4.2) (x^oo) 

provided the integral converges throughout its range for all sufficiently large x. 

The proof of Watson's lemma is a consequence of the fact that the larger part of the area 
under the graph of e~^*, lies over ever smaller intervals of the form [0, e) as x grows. This implies 
that the value of the integral in (4.2) depends mostly on the asymptotic expansion of q{t) when 
t is close to 0, the main term in the formula coming from: 

/ e-^* V aJ'^+^^/'^dt ~ / e-^* V a,t(^+^)/'^dt = V F ( ^ ) 
•^0 s=o -^0 s=o s=o V ^ y 



In our results the functions (/>(t) and q{t) in (4.1) will vary with x. This will require that we 
take care about how things change as x grows. The following result can be found in any book on 
asymptotics, and, is basically a reformulation of Watson's lemma in a form which is convenient 
for subsequent proofs. 

Lemma 4.2. Let {Bn}n>o be a sequence of positive numbers such that Bn — >• +oo. Let 
{4>n{t)}n>o be a sequence of functions, analytic on the unit disk, such that: 

(1) Each 4>n{t) is increasing on [0, 1). 

(2) (pnit) = t -\- 0{t^) as t ^ 0'^ and uniformly for all n. 

Let e > 0. Then, uniformly for all x with Re{x} > e; 



(4.3) 



r ^-B„Mt)t^-idt = s-^r(x) + o(s-^<^^->-i) 

^0 



as n ^ oo. 
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Proof. Since <pn{t) ~ t unifornily for all n ds t — )■ 0, cpni^cx^ > q;/2 for all small enough, ct G (0, 1) 
and for all n. Then, since all are increasing, 



(4.4) Io= f e-'^"'^"(*)t^-^dt < 

J a 



J a 



<ORe(a;}{e ^ }, 



as n — >■ +00. 



(4.5) 



Jo 

^ f g-B„tgB„[t-0„(t)]^i;-l^^ 

Jo 

poo f'OO pOL 

= e-^-H'^-^dt- e-^"H''-Ut+ e--^"*{e^"[*-'^"Wl - 

Jo Ja Jo 



Now, 
(4.6) 



s-^r(x) -h + h 



l-OO 



as n — )• +00. 

By item (2) in the statement of the lemma, f{t) = e^riO{t^) ^ _ q(^j^^ 2 ^ as n — >■ 00, and 

uniformly for all n. Then, there is C > such that f{t) - 1 < e*^-^"*^ - 1 = 0{Bnt^) for t and 

_ 1 

n as before. Taking a = o{Bn ^ ) as n ^ 00, 



(4.7) 



/2 = re--^''*{e^"[*-'^"(*)l - 
Jo 

< S-^<^>-^r(Re{x}) 



as n — >■ 00. 



Finally, take a = (Re{x} + i) '°s-^" in all the above estimates (notice that a is still smaller 
than ^== for large enough n). We get that, as n — >■ 00, 

(4.8) /o,/i,/2«S-^^=^>-i 
Putting all the estimates together: 

(4.9) e-^-'^-Wt^-^di = B--T{x) + ORe{.}{i?n''^"^"'}, 
Jo 



as n — >■ 00. 



□ 
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4.1. The case |r| G [O,logd]. We begin by restricting r to interval [0, logd]. To obtain all 
the subsequent results in the r < case, we just have to notice that Cdix) = Cd{x). A direct 
application of Lemma 4.2 is enough to deal with the present case. We will use a Mellin transform 
representation: 

Proposition 4.3. // Re{x} < 0, then, 

c,ix)=-'^ r £d{-t)t-^-'dt 

Proof. Let B[x, y) = ^-^(x+^y) Beta function. Let {a)j = ^^^^^ be the Pochhammer 

symbol. Then, by well known formulas, 



d 



j=0 
d 



j^f'^^i-iyBij-x^d+i-x) 



j=0 
d 



^ l(^\_iy'[iis^)^(d+l-j + x 



T{d + 1) 

Ti-x)rii + x)±(^'^i-iytMl^ 



j=0 

d 



j=0 



d 



TT -sr-^ (d\(d + X—j 



E 



sin TTX ^ V 7 / V d 

■ . ^ J^d{x) 

sm TTX 



Using the functional equation: 
(4.10) £d{t-') = i-lf +H£d{t) 

We write: 

Cd{-x)= / £d{-t)t--^dt 
Jo 



sm TTX 



J £d{-t)t''-^dt + J £d{-t)t'^~^dt 

C £d{-t)e-'dt + C £d{-t)t(''''^-'dt 

Jo Jo 



About £d{—t) we can say the following: 

(1) £d{—t) is decreasing on [0, 1] for all d. 



□ 
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(2) £di-t) = exp{(ilog(l -t)-{d + l) log(l + t)} = exp{{2d + l)[t + 0{f)]} as t ^ 0+, 
uniformly for all d. 



Then, it is not hard to see that the hypotheses of theorem 4.2 are satisfied, and, for all x G Sd, 
and uniformly for all large enough d, 



(4.11) -^^Cdi-x) 

Sm TTX 



= {2d + l)-^r(a;) + 0[{2d + l)"^"^] + {-lf+\2d + lf--^T{l - x) + 0[{2d + l)^-^! 

4.2. The case |r| — )■ +oo, |t| < d — v^. Next, we will take care of the case |r| G [logti, d- 
The method of steepest descents (two references are |15| and [9]) deals with integrals of the 
form: 

(4.12) I{X) = J e^'^Wg(t)dt, 

where F C C, X is some complex parameter, (/)(t) and q{t) are fuctions analytic in some domain 



containing T. We are interested in what happens to the value of the integral in (4.12) as 
\X\ —7- +00. The idea of the method is to deform T in such a way that I{X) that is amenable 
to the application of Watson's lemma. Roughly stated, the method for achieving this last goes 
as follows: 

(1) Find the saddle points of the function | exp[X(/>(t)]|, that is, where 4>'{t) = 0. Assume that 
a is one of these saddle points, and, the Taylor series of (j){t) at a is (p{a) + Ylm>M Cm{t — 
a)"*, with M > 2 (remember that </>'(a) = 0) and Cm 7^ 0. We know that there are M 
curves Fi, . . . ,Tm which pass through a which satisfy lm.{X(j){t)} = Im{X(^(a)}. These 
last are the curves along which the modulus of e'^'^^^^ varies the fastest. 

(2) Deform F into a new curve F', which passes through the saddle points of | exp[X(^(t)]|. 
F' must also adhere as closely as possible to the "paths of steepest descent" . These last 
are the F^ above for which ^e{XCM} < 0. 

(3) We must check that the size of the integral on sections of F' which are far from the 
saddle points become small as |X| grows. 



Given an integral of the form (4.12 ), it is not at all obvious that the method above can be applied. 
There have been some efforts to make the use of the method easier and more automatic. For 
a useful discussion on this last topic see [13]. In this last the authors advice that instead of 



deforming the curve F in (4.12) onto the paths of steepest descent, we should deform F so that 
it lies on a straight line tangent to the path of steepest descent, that is, integrate over the paths 
of steepest descent of | ex.p[XCM{t — at least when inside the disks of convergence of the 

Taylor series of (j) around the saddle point. Outside of these convergence disks we take paths that 
minimize the absolute value of the integral. In our case this method would work fine for most 
values that t < d could take. Problems do arise when r is close to d. More specifically, when 
d — \t\ = 0(v^). This last is the "coalescing saddle point case" studied by Chester, Friedman 
and Ursell in [8j. In this last case, it is easier to apply slightly modified version of the method 
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outlined by Pagola et al. in |13j . That is, in the vicinity a saddle point a, we integrate instead 
over the path of steepest descent of ew{^['l2m=M C^mit — o)™]}, for some M' > M. 

We remind the reader that t = t^^ = Im{— x^j}. Unless otherwise stated, we assume r > 0. Let 
Qd{t) = log{(^)'^f^}. We look for the zeros of Q^, which are just the zeros of the polynomial 
2dt + irfl — t^). These last are: 



(4.13) 



a± 



V T V T 



d2 
2 



Three things to notice are that 

(1) a± are pure imaginary. 

(2) a+a_ = —1. 

(3) ia_ G (0, (i) and ia„ <C ^• 

Let a = a_. The Taylor series of Qd around a is then Qd{a) + X]m>2 ~ ct)"^) where: 



Cn 



(-1)"*+!^ d (-l)™+iirrf 

(1 + a)™ (1 - a)™ ~^ a™^ 



We can estimate the size of the Cm as follows. Remember that a satisfies 2dt + zt(1 — t^) = 0. 



(4.14) 



Cm. — d 



d 



(-l)™+i(l-a)™ + (l + a)™ (-l)™+iir 



2\m 



(l-a2) 
(-l)"+^(l-a)'" + (l + a) 

2d~ 



+ 



a'" 

.-|^\m+l 



IT 



{-ay 



m— 1 



(-1) 



m+1 



a'" 
I — a 



+ 



1 + a 



a' 



Notice that the bound is uniform in all the a and d, since |a| < 1 and |r| < d. Using the above 
formula, and the fact that a is pure imaginary, we get: 



(4.15) 



Cn 



-H) 



m—l 



\a\ 



where Km is inside the interval (0,2). In particular. 



(4.16) 



d^ 



ASYMPTOTIC DISTRIBUTION OF THE ROOTS OF THE EHRHART POLYNOMIAL OF THE CROSS-POLYTOHB 



Now, for m > 3, 



K — 1 



(4.17) 



> 1 



> 1 



T 

d 
1 

71 



m—l 



(-1) 



m+1 



1 — a 



+ 



1 + a 



2 J V 2 

' 2Re{(i±^^)'"} ; m odd 

2ilm{ (''^^)™} ; "T- even 
m-l „ 



Here we used that < r < d and < ia < 1. Therefore, it is not difficult to see that C3 ^ 



uniformly for all r. At the same time, C2 = o( 



whenever r ~ d as d — )• 00 by (|4.16l). These 



last two facts determine our choice of integration contours. 



Since £d satisfies the functional equation (4.10), and is real on the real line, we can restrict 



ourselves to working inside D. This last will become clearer as we go through the calculations. 
We first define the restriction of our integration path to D: 



A. 



The reason for the subscript d is that the path will actually vary with d. First, we de- 




FiGURE 3. Integration path r^; restricted to D in solid black. The actual steepest 
descent path is in gray. The cubic through the saddle point a and tangent to the 
steepest descent path in dashes. 



fine the section of path which passes through the saddle point a and is tangent to path of 
steepest descent that goes through a. We choose the path of integration in such a way that 
lm{^2^(t - af + ^{t - af] = and Re{^(t - a)^ + ^(t - a)^} < 0. Now, from \A.IQ\ and 
( |4.15 ) we have: 



C2 
C3 



I a 



1^ 
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We finally define: 



Aa(r) = a + zgr 



T + 1 
3T + 1 



+ ="+2H^^ 



r + 1 
3r + 1 



We will determine the domain of Aq, below. This last will be a subset of (— |,oo), below. We 
thus have: 



2 



(A«(r)-a)2 + ^(A,(r)-a)3 



C2 

2 



3, |-f^2\ ^2 



8 



T + 1 
3T + 1 



T + 1 
3T + 1 

\ 2 



+ i + 



C3 /3, ,i^2 



2' 'K3 



T + 1 

3r + i 



= _9^^j.2(22:+i)VrTT 

4^i^| (3r + l)i 

From the above equation define: 



T+1 
3T + 1 



+ ^ 



(4.18) 



F{T) := 2T 



2(2r + i)^VrTi 

(3r + l)i 



F(T) is positive on (— g,oo), decreasing on (— 1,0) and increasing on [0,oo). To see this last 
just write: 



F'(T) 

F'iT) = F{T)^ = F{T) 



2 4 I I 

r'^2r + i'^r + i 3r + i 



and, use the fact that 2r + 1, T + 1 > 3T + 1 > when T < 0, while < 2r + l,r+ 1 < 3T + 1 
when T > 0. 

Write (^[A„(r)] = -|r§F(T) - G(r). We want Rej It^F{T) + G(r)} to grow like F{T) 
in the domain of definition of Aq,(T). Assuming we have chosen the domain of AQ,(r), we have: 



/ £d{t)t^-^dt = 

J Ac 

-L 



exp{Q(i)}- -dt 



3| 1^2^ , , 



q TV'S 

-1 / exp^--r-^F(r)-G(T) 



2 A3 



i + |M||c/(r) 



(T-l 



-U'{T)dT 
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Now, from (|4.15l), 



\m—l 



G{T) = Y.^ [A„(T)] 

m>4 

= E -(-) 

m>4 



-Kr, 



m>4 



3, |i^2^ 



T + 1 

3r + 1 



+ i 



r + 1 

3T + 1 



Let U{T) = T [\l Then, \U{T)\ = Then, assuming \U{T)\ is smah 



enough, for example, if < 
(4.19) 



G(r) 



m V2^3 



U{TY 



2kJ 



U{T)^{1 + 0[K2U{T)]} 



Then, for smaU enough T, 

(4.20) exp{-G(r)} = 1 + 0{tKIu{T)^} 

Using that ia G (0, 1), and a varies in a compact subset of (0, 1), it foUows that, for smah enough 
T, it is not hard to see that: 



(4.21) 



1 + 



3 H 
2 a K3 



U{T) 



cr-l 



1 + 0{K2U{T)} 



1 _ ^\^K2JJ^r|.^ 
^ 2 l-a Ko,'-^ ' 

It follows that for small enough T, we can write as: 

(4.22) Za{x) = l^exp |-^r^F(T)| U' {T) {l + O [tKI\U{T)\' + K2\U{T)\\ } dT 



We choose J so (|4.19|), (|4.20|) and (|4.2lD are satisfied; 

1 



(4.23) 



J 



1 



1 + ^-4/7^-2 P^2/7_^, 



1 



Then, for T G J, |C/(r)| = 0{t^^^'^ K^'^) = oir-'^/^K^^). Now we can write ( |4.22| ) as: 
(4.24) X,(x) = |^exp|-^T||F(T)|c/'(r)dr{l + 0[r-i/28]} 

This last estimate is sufficient for the purpose of finding the zeros of Cd- The estimate we obtain 
below is cleaner looking, and also suffices for the zero search. Now, we estimate the tails. F(T) 
is increasing on [0, 00), and, by (4.18), F{T) S> for all T > 0. Also, 



(4.25) 



U'{T) = i + T 



T + 1 r 1 1/2 3/2 



+ 



3T + 1 T T + 1 3r + 1 



o(i), 
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as T — )• oo. We then have: 



(4.26) 



9 Kl 

^2 
'■3 



exp<J-^rpF(r)|>[/'(r)dr« 



POO 

exp {-ri^lCT} T^'^'^dT 



3 Jr-6/7i^-3 



as d — )• oo, for some positive constant C that does not depend on d. Since we assumed < 
d — T = il(v^),we get: 

(4.27) r exp\--T^F{T)] U' {T)dT « exp{-Cri/^} 



From (4.25) and (4.18), we get: 
(4.28) 



U'{T) < (ST + l)-3/2 ^ 



for all T G (-^,0]. It follows, then, 



(4.29) 



1/3 



9 Kl 



-1/3 



1+^-4/7^-2 



9 '■^ 
2 



exp ^-tKId{3T + 1)~3/2| (32^ _^ ly^/'^dT 
exp{-TKlDT} T-^/^dT 

< exp{-L»r^/^} 



[3(l+^-4/7^-2)]3/2 



as (i — )• OO, for some positive constant D that does not depend on d. Putting together ( |4.24[ ), 



( 4.27| ), and (|£29j), we finally get: 

/oo r Q "I 



(4.30) 



as d — >• oo. The above equation is the main term of the asymptotic formula we are looking for. 
Remember that we started by finding a saddle point of (j)d, a G D. Then, we integrated £d{t)t^~^ 
over small curve segment that is tangent to the path of steepest descent through a. This last 



was Aa{T) = Q + llall^T' ^r+i + ^j- For our choice of domain J it isn't hard to see that, 
for T € JnM>o, Ao,{T) = a + Oir^"^/"^) will be contained in the fourth quadrant of the complex 
plane, while for T £ J f] M>o, Aa(T) will be contained in the third quadrant of the complex 
plane. This last will give the main term that comes from the vicinity of the saddle point a. In 
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reality we want to integrate over — Aq,. So, the contribution to our asymptotic formula coming 
from the vicinity of a is, as d — )• oo: 
(4.31) 

o roo ( n r^3 'I , ^ 

-1/28 



£:,(t)t-<*-idt = -^^^^:,(a)a^y""^exp|-^T^F(T)|c/'(r)^^ 



Inside the unit disk, the rest of our integration will be done over line segments. We will choose 
these line segments in such a way that the value of the integral over them will be small. We want 
to make sure that over the line segments the modulus of exp Q(t) decreases, as t moves away 
from a, and is bounded above by the modulus of exp Q(a). One way of doing this is by noticing 
that if t G D has Re{t} < 0, then, the modulus of SaiTt) decreases, while that of (TtY'^ remains 
constant, as T goes from 1 to l/\t\. In other words the modulus of £d{Tt) decreases as as T 
goes from 1 to l/\t\. A similar analysis shows that if Re{t} > 0, then, the modulus of £d(Tt) 



increases as as T goes from to 1. We can finally write down all of n D. Write J (see (4.23) 
)as [a_,a+]. Let /_ be the line parametrized as /_(T) = TAQ,(a_), where T G [1, |AQ(a_)|~^]. 
Let /+ be the line parametrized as 1+{T) = TAQ,(a+), where T G [— | Aa(a+)|~^, | AQ,(a+)|^-'^]. 
Then, our integration contour inside D will be: 

Arf = /+ - A„ + l_ 

One should note that actually touches the origin, which in principle is not allowed, since F 



should be as in (1.3). To see that this last does not pose a problem, just notice that integral of 
£d{'t)t^~^ on a straight line through the origin goes to as the length of the line goes to 0. To 
bound the values of the integrals over and I- we first note that: 

(max I exp Q|)(min|l-t|) / |tr^-V*l < 

max I exp Q\ 

Ji± '± '± Ji± '± 

Here we used the fact that the distance from l± to 1 is bounded away from zero and that the 



integrals J^^ Itl"^ are uniformly bounded above. Now, by (4.20) and (4.28), 



exp[Q(a±)] = exp[Q(a)]exp{(/)[A„(a±)]} 

9 KS 



( 9 K 

exp[Q(Q)] exp |--r-^F(a±) - G{a±) 
exp[Q(a)] exp |-^r^F(a±)| {1 + 0(rK|t/(a±)4)} 



9 Kl 

= exp[Q(a)] exp{-Z)r^/^}{l + 0{t-^/^^)}, 
where D is some positive constant. Then, 

/ £d{t)t''-^dt < |A„(a±)|exp[Q(a)]exp{-L>ri/7} < £d{a)a'^+^ exp{-DT^/^} 
Ji+ 
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We then have: 



/ £d{t)e-^dt =( [ + [ + [) SSW'^dt 
JAa \J-Aa Jl- Jl+J 



(4-32) = _i^:p^,(«)a-I(a){l + 0[K^' eM-Dr'/')]} 

= ~i~£dia)a-Aa){l + 0[exp(-£)ri/^)]} 

In the last Hne above we used that a £ [e, 1 — e], and that K2 = -^1 — T^/d"^ 3> d~^^^'^ since 
< r < d — v^. Finally, we define the section of integration contour that lies outside B. For 
any given curve in 7 : [0, 1] ^ C define 7"^^) = [7(^)]~^ = j{T) and -j{T) = 7(1 - T). 
Finally, define F^ n (C \ B) := - {TJnWj'^. We then have: 



/ Sd{t)e''-'dt= [ ^d[-(rdnB)-i][-(rdnB)-i]-'^-H-i)[-(r,nD)-2][-(rrfnD)r 

^-(FdnD)-! J 

= {-if j £d[-{^d n B)][-(rd n B)]-^[-(rrf n d)]' 

= f £d{t)t(^-^)-^dt 

Note that Im{l — x} = r, so a and cf) do not change. We can then use all the above estimates 
to obtain, as d — >■ +00: 

(4.33) / £d{t)e-^dt = [ + [ 



/rdn(C\D) 



£d{t)e-^dt = 7{d, x) + {-lY+^7{d, 1 - x), 



where, 

(4.34) F{d,x) = -Z2^^<i(«)«^Xrf(a;){l + 0[K2"'exp(-Z^Ti/7)]}, 

where, 

9 



i{x) = |"^exp|-^r^F(r)|c/'(T)dr 



2d I 

■ -1/3 L ^ ^3 
Note that -iZK2/2K^ = 3C2/2C3. 

5. Some commentary 

The asymptotics of Cd{x) can be written down for all x with Rc{x} = —1/2. Remember 
that we worked only enough of the asymptotics of Cd to find its roots. The methods used here 
seem to be perfectly applicable to the study of other polynomials and functions that have a 
nice enough generating series. For example, if we take some sequence of polynomials {Vd}d>0: 
where d is the degree of Vd-, such that the generating function Ylim>o'^d{rn)t"^ has all of its 
zeros on the unit circle, then, the methods used here can be applied under the appropriate, 
not very strict conditions. These last topics will be subject of follow-up work. Much of the 
work presented is inspired on the methods used in the theory of orthogonal polynomials and 
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the theory of analytic combinatorics (for a very thorough introduction to the field of analytic 
combinatorics see Flajolet's and Sedgewick's book [lO]). The polynomials Cd form a family of 
orthogonal polynomials which happen to be Mellin transforms of Laguerre polynomials (see [6] ) . 
In general one can construct, through the use of recursions, families of orthogonal polynomials, 
which have generating series with nice properties, such as having generating functions with zeros 
on the unit circle. The methods used in this paper can be applied to this last as well. On the 
combinatorial geometrical side of things there might be interesting applications to knowledge 
about the roots of Ehrhart polynomials; see Henk's et al. work in [TT]. 
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